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Abstract 



We present a numerical method for generating vortex rings in Bose- 
Einstein condensates confined in axially symmetric traps. The vortex ring 
is generated using the line-source approximation for the vorticity, i.e., the 
rotational of the superfluid velocity field is different from zero only on a cir- 
cumference of given radius located on a plane perpendicular to the symmetry 
axis and coaxial with it. The particle density is obtained by solving a modi- 
fied Gross-Pitaevskii equation that incorporates the effect of the velocity field. 
We discuss the appearance of density profiles, the vortex core structure and 
the vortex nucleation energy, i.e., the energy difference between vortical and 
ground-state configurations. This is used to present a qualitative description 
of the vortex dynamics. 
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I. INTRODUCTION 



Since 1999, when vortex lines in a trapped Bose-Einstein condensate (BEC) were first 
experimentally obtained their study has received great experimental and theoretical 
interest as it constitutes a clear signature of superfluidity effects in these confined systems. 
Some remarkable experimental achievements are, among others, the study of the dynamics 
of single vortex lines 0] and the formation of small and large [Q vortex arrays. A 
review of the research done in this field is presented in Ref. [Q. It is worth to note that 
the experimental situation in BEC regarding the formation and detection of vortices is at 
variance with that in helium II droplets. For these drops, a paradigm of superfluid finite 
systems, the problem of nucleating vortices and their stability has only been addressed very 
recently from the theoretical point of view (see Refs. and their experimental detection 

is still an open question. 

Vortex rings are vortices whose core is a closed loop with quantized circulation around it 
||. They are complex topological structures that have attracted and will continue to attract 
some experimental and theoretical interest. Based on numerical simulations, different meth- 
ods have been proposed to generate vortex rings in BEC. As in bulk liquid helium, vortex 
rings may be produced introducing an impurity in the condensate with a definite velocity, 
whose displacement causes a vortex ring ||1U|| . Another mechanism [II] consists in using dy- 
namical instabilities in the condensate to cause dark solitons to decay into vortex rings. A 
well controlled method to produce vortex rings by electromagnetically induced atomic tran- 
sitions in two-component condensates has also been put forward ||12|| . The method proposed 
in Ref. has been successfully applied in Ref. fl3| to generate vortex rings experimentally. 

Rather than proposing a method that could be implemented to create a vortex ring 
experimentally, our aim here is to set up a numerical method simple yet accurate enough to 
generate quantized vortex rings with a definite radius R in one-component condensates. We 
restrict our study to vortex states such that the divergence of their velocity is vanishingly 
small. This assumption yields analytical expressions for the velocity field, and provides 
a fair approximation to the superfluid flow around the vortex core for rings in the bulk 
of large condensates (Thomas- Fermi limit). We have considered large condensates at zero 
temperature and therefore, dissipation has not been taken into account. They are axially 
symmetric about the z axis, and have the z = plane as symmetry plane, and may host a 
vortex ring of radius R, coaxial with the symmetry axis of the trap, and placed on a plane 
at a distance Z > from the symmetry plane. Although generalization to situations with 
more than one such vortex rings is straightforward, we have not considered this possibility. 

This paper is organized as follows. In Sect. II we present the method used to generate 
vortices. The way to obtain the velocity field of a vortex ring is described in Sect. Ill, 
and the explicit expressions of the velocity components are given in an Appendix. Sect. 
IV is devoted to the analysis of the particle density profiles and vortex nucleation energies, 
which allows a qualitative description of their dynamics. We also present results obtained 
from a completely different, more involved method we have set up to generate ring vortices 
that permits to test the approximation of zero divergence of the velocity. Finally, a brief 
summary of the results is presented in Sect. V. 
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II. ENERGY FUNCTIONAL 



We consider a weakly interacting Bose-condensed gas confined in a harmonic trap V ex t(r) 
at zero temperature. In the Gross-Pitaevskii (GP) theory, the ground state (g.s.) energy of 
the condensate is given by the functional |14| 



dr 



'*|V»| 
2m 



+ V C 



cxt I 



l^rl 2 _|_ £|^r| 4 



(1) 



where ^(r) is the condensate wave function. The first term in Eq. (ffl) is the kinetic energy of 
the condensate, the second term is the harmonic oscillator energy arising from the trapping 
potential, and the third term is the mean-field interaction energy. The coupling constant 
is g = Airh 2 a/m, where a is the s-wave scattering length, and m is the atomic mass. The 
number of atoms in the condensate is /<ir|\I/| 2 = N. The g.s. wave function is determined 
by solving the GP equation obtained minimizing the energy functional. 



The wave function \I/ can be written in terms of the particle density p(r) 
phase S(r) as 



l#( 



(r) = Jp(r) exp[iS(r)] 



and 



(2) 



and the superfluid velocity is given by v = (h/m)'VS. In this work we will use the equivalent 
quantum hydrodynamic description of the condensate in terms of the density and the super- 
fluid velocity |14[], since it allows a straightforward generalization of the energy functional 
Eq. (0) to include vortex states. Using Eq. (|2]) it follows that 



E[^] = E [p] + EUP^ 
The first term is only density- dependent: 

n 



E [p: 



dr 
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(3) 



(4) 



The first term in E [p\ is the quantum kinetic energy. The second term of the energy 
functional Eq. (|3|) corresponds to the kinetic energy associated to the flow velocity v, and 
is given by 



-m 



dr p(r) 



(5) 



The g.s. wave function in the absence of vortices has a spatially constant phase and therefore 
zero velocity. It can be obtained by minimizing Eo[p\. We consider condensates with 
positive scattering length and axially symmetric traps V r ext (r) = m[u) 2 ± (x 2 + y 2 ) + u 2 z 2 ]/2 
[= m(uj_r 2 +uj 2 z 2 )/2 in cylindrical (r, z) coordinates], with different values of the asymmetry 
parameter A = uj z /uj±. The trap harmonic frequency uj± provides a length scale for the 
system, a± = (h/muj^) 1 ^ 2 . We will use a±, Hu±, and N/a]_ as units of length, energy, and 
density, respectively. 

If we consider large condensates in which the Thomas- Fermi (TF) approximation holds 
14j, the quantum kinetic energy can be neglected compared to the interaction energy in 
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Eq. (£|), and the g.s. density of the condensate in the absence of vortices is given by 
Po(r, z) = n(l — r 2 / R\ ¥ — z 2 /Z^ F )/g in the region where this expression is positive, and 
zero elsewhere. The TF extents of the condensate in the radial and axial directions are 
i?TF — (2p/mtu 2 J 1 / 2 and Ztf = Rtf/^- The chemical potential p is fixed by normalization 
p = hu±(l5\aN/ a±) 2 / 5 /2. We recall that the validity of the Thomas-Fermi approximation 
is guaranteed if Na/a± 3> 1. 

We turn now our attention to the case of condensates with vortex states characterized 
by a given irrotational velocity field associated with a non-vanishing quantized circulation. 
The total energy of the system is given by the energy functional Eq. (^). If the velocity field 
is known, for a given number of particles the density profile can be obtained minimizing 
Eq. (H). This yields the equation 

/ ft 2 v 2 1 \ 

+ Kxt(r) +g |^(r)| 2 +-mv 2 U(r) =/i^(r) , (6) 



V 2m v ' ° 1 r w 1 2 

where ip = yfp is the modulus of the complex wave function Eq. (0). This is the GP equation 
expressed in terms of the hydrodynamic variables. 

In the presence of a quantized vortex, the density of the system drops to zero at its core, 
whose size is characterized by a healing length £. For very large condensates the healing 
length can be approximated by £ = (Sixp^d)^/ 2 , where po is the density of the condensate 
before creating the vortex. For a centered vortex line in the TF approximation p = p/g, 
and the corresponding healing length £ is 



Rtf V Rtf 



(7) 



with £o <C a± <C Rtf- In this approximation a local healing length can be defined |Tj| as 

£(r, z) = £o/ yl — (^/-Rtf) 2 — (-z/^tf) 2 - Note that the size of the core is larger for a vortex 
in the low-density region. 



III. SUPERFLUID VELOCITY FIELD 

The velocity field around a straight vortex line has an analytical expression when the 
vortex is along the symmetry axis 0,|T^]. Approximate analytical expressions can be found 
in the case of vortex lines off the symmetry axis in large and very elongated, quasi-two 



dimensional condensates (see Ref. |17] and references therein). However, vortex lines gener- 
ally bend in three-dimensional condensates P,l8| which renders impracticable an analytical 
treatment of the velocity field. 

In the case of quantized vortex rings we proceed as Schwarz and Jang have done in the 
case of helium II |19| (see also Ref. |p0[). For a vortex ring characterized by the values of 



(R, Z) already defined, circulation number n = 1, 2, . . . and quantum circulation ko = n h/m } 
we write the vorticity in the line-source approximation: 

w = ko8(r-R)6(z-Z)4> , (8) 

where (r, z, <p) are the cylindrical coordinates, and <p is the unit vector in the azimuthal 
direction <ft = (— sin 0, cos0, 0). The superfluid velocity field that arises from this distributed 
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vorticity fulfills w = Vxv. Hence, the velocity field around the vortex is irrotational except 
on the vorticity line, where the density of the condensate is zero. If V • v = to a good 
approximation a velocity vector potential A(r) can be introduced such that v = V x A. 



If the vorticity is specified, A(r) is determined by the equation 

uj = V x (V x A) , (9) 



whose integral solution is |1| A(r) = A (r, z)4>, with 



At \ k ° O COS0U/)' 

A o(r,z) = —R (10) 
47r Jo ^r 2 + R 2 -2r J Rcos0' + (z - Z) 2 

The radial and ^-components of the velocity are obtained as 

v dA ° 

dz 

(11) 

1 9 / AS 

r or 

whereas the azimuthal component of the velocity field around the vortex ring Eq. (§) is 
zero. We give in the Appendix the general expressions of Aq, v r and v z written in terms of 
hypergeometric functions |22] . 



For a confined system, the existence of a boundary and the fact that the density is 



inhomogeneous may have an effect on the actual velocity field p3fl . In our case, since the 
condensate extends up to 'infinite distances', there is no need to introduce any image vortex 
to ensure that there is no particle flow across the boundary due to the superfluid motion. 



Moreover, it has been shown [15j that in the TF limit the corrections to the velocity due 
to density inhomogeneities can be safely neglected. Thus, we approximate the velocity 
field by Eqs. ( |T0|) and fllTD (see Appendix), but have restricted ourselves to study vortex 
ring configurations (R, Z) whose vorticity line is inside the domain where the density is 
positive in the TF approximation, i.e., only vortex rings (R, Z) satisfying the condition 
(R/Rtf) 2 + {Z/Ztf) 2 < 1 are considered. We call the (R, Z) line defined by the condition 
(R/Rtf) 2 + (Z/Z TF ) 2 = 1 the TF boundary. Once the the superfluid velocity field has been 
fixed, the density profile of the condensate is obtained solving Eq. (|6|). 



IV. RESULTS 

We present numerical results for large condensates in the Thomas-Fermi limit hosting 
a quantized vortex ring with circulation number equal to one (nucleating vortex rings with 
n > 1 is energetically less favorable P,p4l). For a singly quantized vortex ring configuration 
(R, Z) with vorticity given by Eq. ® and circulation number n = 1, we have computed 
the velocity field, Eq. ([11]), and have obtained the density profile of the confined condensate 
by solving Eq. (||) in cilyndrical coordinates using the imaginary time method [25]. Note 



that although we are in the TF limit, we do solve the complete GP equation to obtain the 
density profiles. 
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Figure [I] shows density profiles in the z = plane as function of r for the experimental 



parameters of Ref. [13], that is, N = 3 x 10 5 atoms of 87 Rb (scattering length a = 5.82 x 
10 -9 cm) in a spherical trap with uj±/2tt = 7.8 Hz. We have plotted two configurations 
with a vortex ring located in Z = having a radius R = 2a± (dotted line), and R = Aa± 
(dot-dashed line), respectively. The density profile of the condensate without vortex is also 
shown (solid line). 

The density is zero on the vorticity line, as seen in the density profiles. As expected, the 
size of the core is of the order of the healing length. Indeed, for this condensate one has 
R TF ~ 5.8 a± and £ ~ 0.17aj_. Assuming that the vortex diameter is twice the local healing 
length, and using the TF expressions we get as core diameters the values 2x£(2, 0) = 0.37 a±, 
and 2 x ((4,0) = 0.47 a±. One can see from Fig. |l| that the core sizes are in agreement with 
these estimates. When the radius of the vortex ring increases, the core lies at a lower density 
region and therefore its size increases. 

To study the vortex ring energetics we have chosen a larger condensate made of N = 10 6 
atoms of 87 Rb confined in an axially symmetric trap with axial frequency uj z /2h = 220 
Hz and three different geometries, namely spherical, disk-shaped and cigar-shaped, with 
A = 1, and 0.2, respectively. For this condensate, we show in Fig. |2] several equidensity 
lines (arbitrary values) in the y = plane. The top panel corresponds to the cigar-shaped 
trap (A = 0.2, i? TF = 8.7 a±, Z TF = 43.3 a^), the middle panel to the spherically symmetric 
trap (A = 1, -Rtf = Z TF = 10.2 a±), and the bottom panel to the disk-shaped trap (A = V», 
Rtf = 45.5 a_i_, Ztf = H-3 a±). The x and z axes are in units of a±. For all geometries, the 
vortex ring configuration is (R = 3.1 a±, Z = 0). The intersection between the vorticity line 
and the y = plane appears as two dark dots, indicating the steep density depression around 
the vortex core. The presence of the vorticity causes a drastic distortion of the density with 
respect to the g.s. profile. Yet, for large condensates it is a rather local effect, as can be 
seen from Fig. |2|, and also from Fig. [I] for the smaller condensate. 

The nucleation of a vortex has an energy cost, since the energy of a condensate with 
a vortex is always larger than the energy of the condensate without it, Eq$. The vortex 
nucleation energy, E — Egs, corresponding to the condensate of the spherical trap in Fig. 
is plotted in Fig. |3| as a function of Z for different values of R. Fixed R, the nucleation 
energy is maximum when the ring is in the z = plane, and it decreases as Z increases, 
since displacing the ring outwards implies that the superfluid flow affects less atoms in the 
condensate. Note that to have the vorticity line within the TF boundary, the .R-curves end 



at different values of Z/Z TF given by y 1 — (R/Rtf) 2 - The nucleation energy does not vanish 
when the ring is located on the TF boundary as there still exists a superfluid velocity field 
which produces some effect. Since we treat the superfluid velocity as an external field, the 
nucleation energy of very superficial vortex configurations may be somewhat overestimated. 
This drawback is apparent beyond the TF boundary, and shows up as a long-tailed nucleation 
energy as a function of either Z/Z TF or R/R TF , although it goes eventually to zero. 

The nucleation energy is plotted in Fig. |as a function of R for different Z values. The 
smaller ring we have considered has Rj R^f ~ 0.1. The middle panel corresponds to the same 
spherical condensate as in Fig. |^. Besides, the top panel corresponds to the cigar-shaped 
trap, and the bottom panel to the disk-shaped trap. The top curve in all panels corresponds 
to vortex rings with Z = 0, and for a given trap geometry (fixed panel), the higher the curve, 
the lower the Z. As in Fig. || we have only described vortex ring configurations inside the 
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TF boundary. For this reason, fixed a Z value the corresponding nucleation energy curve 
stops at the value R/Rtf '■ 



>i-(z/z TF y. 

Figure |] shows that for each asymmetry parameter A, given a Z there exists a vortex 
configuration corresponding to a certain radius R that maximizes the nucleation energy. 
This radius increases as Z decreases. In the (R, Z) two-dimensional configuration space, the 
nucleation energy E(R, Z) — E GS has only one maximum at (R eq / R^f — 0.6, Z = 0). It is 



in agreement with the results obtained in Ref. p6| for a spherically symmetric trap using a 
simplified model of a vortex ring. 

This maximum corresponds to the only stationary, although unstable, vortex ring con- 
figuration (equivalent to a vortex line nucleated along the z-axis), whose location is almost 
independent of A but its value depends on it, decreasing as A does. This means that the 
energy to nucleate a vortex ring characterized by the values of (R, Z) is lower the more 
elongated the trap is. 

For each asymmetry parameter, the 'iso- nucleation energy' lines in the (R, Z) configura- 
tion space corresponding to a value not too far from the maximum value are closed curves 
around the maximum, and the allowed ring radii are in the range defined by the intersection 
of a straight line representing the given nucleation energy and the Z = curve in Fig. |]. If 
we relax the condition that the vorticity line must be inside the TF boundary, this should 
be always like this ||26|| . Otherwise, the isoenergy lines eventually do not close. Yet, the 
allowed ring radii are determined in a similar way. 

We have considered the vortex ring as if it were a static object. However, in a trapped 
condensate a vortex ring will move with a velocity that results from the interplay between 
the effect caused by the inhomogeneity of the condensate and self-induced effects arising 
from its own local curvature |||9],^6|. A qualitative analysis of the dynamics of coaxial 
vortex rings can be carried out from Figs. § and f|, as proposed in Ref. , and tested in 
a real time dynamics for two-component condensates [ff^LBBfl . It can be seen that a coaxial 
vortex ring oscillates inside the condensate, moving towards the surface and instead of being 
annihilated at the boundary, the ring generates a background flow in the rest of the cloud 
that draws it back towards the center. The ring dynamics in a trapped condensate results 
in an oscillation along the symmetry axis of the trap. The radius of the vortex ring grows 
or shrinks depending on its Z-position, in such a way that the ring motion almost follows a 
trajectory with constant energy in the (R, Z) configuration space. 

Finally, we have used a completely different method, numerically more involved, to gen- 
erate quantized vortex rings. It does not make any assumption on the superfluid velocity 
field, nor on the radius or position of the vortex ring either, and can be used to test the 
results we have obtained based on the validity of the approximation V ■ v = 0. 

We have proceeded by analogy with the case of a vortex line, and have added a constraint 
to the energy functional that forces to form a toroidal-like hole in the condensate with a 
flux of atoms around its core with zero azimuthal velocity component. We have used as 
constraining operator the angular momentum about the axis, i.e., 0-L = — ih(rd z —zd r ) 



27 



Introducing the functional 



dr 



2m 



\VW + V ext (r)\V 



+ g -m 

2 1 1 



+ fy / dr^L^ 



(12) 
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where the Lagrange multiplier can be understood as the local angular velocity about 
the azimuthal direction, we have obtained the associated GP equation. For a given Q^, 
we have solved the partial differential equations obeyed by the real and imaginary part of 
the condensate wave function, which are coupled by the presence of the constraint. These 
equations have been discretized on a (r, z) mesh using seven-point formulas to represent the 
differential operators, and have been solved using again the imaginary time method | f25H . 

Taking as an example Q$ = 0.2 uj±, we have obtained a coaxial vortex ring configuration 
with its vorticity located on the z = plane whose equidensity lines cannot be distinguished 
from these represented in the middle panel of Fig. |2| at the scale of this figure (whose 
parameters were chosen before to allow us to make this statement). 

Figure presents a comparison between the two methods. It shows the particle density 
profile as a function of r at z = for the spherical condensate of Fig. |^. The solid line is 
the configuration that minimizes Eq. fljjD - The density goes to zero in R — 3.1 a±, and the 
vortex ring has an energy equal to 37.164 hu±. The dashed line is the configuration obtained 
by the method of Sect. 2 with the ring vorticity Eq. (§) placed in (R = 3.1 a±, Z = 0) which 
has an energy of 37.167 tku±. Thus, the energy and density profile of both configurations are 
almost identical, and the velocity fields are also in good agreement. 

We have also checked that minimizing the functional Eq. ( |T2"| ) yields quantized ring 
vortices with n — 1. To this end, we have proceeded to a direct numerical integration of the 
circulation § v • dl using several arbitrary closed paths in the (r, z) 'plane'. We have found 
that the circulation for paths enclosing the vortex core is h/m with a good accuracy, and 
is zero otherwise. The velocity field has been calculated from the wave function recalling 
that the current field is j(r) = pv = {Ti/2mi)[^!*V^ - (V\l>*)\&]. The current field j(r) 
in the y = plane is plotted in arbitrary units in Fig. |6|. The axes are in units of a±. It 
can be clearly distinguished in this figure the position of the vortex ring core at Z = and 
R ~ 3.1 a±. 

V. SUMMARY 

We have developed a method to study vortex rings hosted in confined condensates in 
the TF limit. It provides analytical expressions for the superfluid velocity field that satisfy 
the irrotational condition and the quantization of the circulation of the superfluid flow. The 
velocity around the vortex core is introduced as an external field in the Gross-Pitaevskii 
energy functional, which is minimized to obtain the wave function of the condensate. 

Using this method, we have computed the density profile and nucleation energy of vortex 
rings in one-component condensates. We have found that the size of the vortex ring core is 
of the order of the healing length as in the case of vortex lines. The presence of the vortex 
causes a sizeable distortion of the density, as atoms are pushed off the core, but the effect is 
rather local. 

The analysis of the dependence of the nucleation energy on the position and radius of 
the ring allows one to conclude that a Hamiltonian dynamics will lead to oscillations of the 
vortex ring along the symmetry axis of the condensate, increasing and decreasing its radius 
in accordance with previous predictions. A dissipative dynamics would cause the vortex to 
decay, for example disappearing across the border of the condensate, as occurs for vortex 
lines 



S 



The method provides a handleable way to generate vortex rings in confined condensates 
that can be used as a starting point to study the dynamics. A detailed calculation of the 
vortex ring dynamics in a confined condensate, either Hamiltonian or dissipative, is beyond 
the scope of the present work and will be developed elsewhere. 
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APPENDIX: 



In this Appendix we give the expressions and method we have used to obtain the velocity 
field. From Eqs. ( [TTj| - |TT|) one has: 



v r (r, z) 



k R(z- Z) 



71 



[{r + R) 2 + {z- Zff 2 Jo {1-8 cos 2 0) 3 / 2 



*/2 2 cos 2 . 



(Al) 



where we have defined s = ArR/ [(r + R) 2 + (z — Z) 2 ]. The integral Eq. ( |AT| ) is written in 
terms of hypergeometric functions |23 yielding: 



v r (r,z) 



k R(z- Z) 



2 [(r + R) 2 + {z-Z) 



213/2 



[F(3/2,3/2;2;s)-F(l/2,3/2;l;s) 



(A2) 



Similarly, 

v z (r,z) -- 



k^R 



+ 



2r[(r + R) 2 + (z - Z) 2 ^ 1 ' 
k R 



m [F(3/2,l/2;2;s)-F(l/2,l/2;l;s)] 



2 [(r + R) 2 + (z- Zff 12 
3R 

+ — F(5/2,3/2;3;s)] 



[(r + R) F(l/2, 3/2; 1; s) - (r + 2i2) F(3/2, 3/2; 2; s) 



(A3) 



and 



A (r,z 



k n R 



2 [(r + R) 2 + (z - Z) 



172 [^(3/2, 1/2; 2; s) - F(l/2, 1/2; 1; s)] . (A4) 



21 V 



The hypergeometric functions entering Eqs. {~K2-~KA) can be written in terms of the complete 
elliptic integrals E and K [^] as follows: 



F(l/2,1/2;1; S ) = -K(v^) 

7T 



(A5) 



F(l/2, 3/2; 1; a) = F(3/2, 1/2; 1; s) = - E< V S ' 



7T 1 — S 



(A6) 
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F(3/2, 1/2; 2; s) = F(l/2, 3/2; 2; s) = — [K(^) - E(^)] 

7TS 



(A7) 



F(3/2,3/2;2;s) 



71S 



1-s 



K(V5) 



(A8) 



F(5/2,3/2;3;s) 



16 

37TS 5 



2 - s 



1 - s 



E(v^) - 2K(v^) 



(A9) 



To evaluate E and K we have used polynomial approximations |2Tjl . 
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FIGURES 




r (a.) 



FIG. 1. Density profile as a function of r at z = for a condensate in a spherical trap with 
the experimental parameters of Ref. (aN/a± = 440) hosting a singly quantized vortex ring at 
Z = 0. The dotted line corresponds to a vortex ring with R = 2a±, and the dot-dashed line to a 
vortex ring with R = 4aj_. The solid line is the g.s. density profile in the absence of vortices. The 
densities have been normalized to one. 
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FIG. 2. Vortex ring state equidensity lines, in arbitrary units, plotted in the y = 
plane (axes are in units of a±). The ring is at (R = 3.1a±,Z = 0) in a conden- 
sate with N = 10 6 atoms of 87 Rb, confined in an axially symmetric trap with axial fre- 
quency lu z /2it = 220 Hz and different trap geometries. The top panel corresponds to a 
cigar-shaped trap (A = 0.2, i?xF = 8.7aj_, Ztf = 43.3aj_), the middle panel to a spheri- 
cally symmetric trap (A = 1, Rtf = 10.2 aj_), and the bottom panel to a disk-shaped trap 
(A = y/8, R TF = 45.5 a ± , Z TF = 11.3 a±). 
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FIG. 3. Nucleation energy of a vortex ring (R, Z) as a function of Z for different values of the 
radius. The labels of the curves are in units of a±. The condensate has N = 10 6 atoms of 87 Rb 
and it is confined in an spherically symmetric trap with frequency u> z /2ir = 220 Hz. The TF radius 
is i?TF = -^tf = 10.2 a±. 
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FIG. 4. Nucleation energy of a vortex ring (R, Z) as a function of R for different values of 
Z. The condensate and trap geometries are these of Fig. ||. The top panel corresponds to the 
cigar-shaped trap, the middle panel to the spherically symmetric trap, and the bottom panel to 
the disk-shaped trap. From top to bottom, and in units of a±, the curves correspond to Z = 
to 40 in Z-steps of 5 (top panel); Z = to 10 in Z-steps of 1 (middle panel); Z = 0, 1,2, and 3 
(bottom panel). 
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FIG. 5. Density profile as function of r at z = for the N = 10 condensate in a spherically 
symmetric trap with a singly quantized vortex ring. The solid line corresponds to the density 
profile that minimizes the functional Eq. ( |l2| ) with constraint Q$ = 0.2 u>±. The dashed line is 
the density obtained imposing a ring vorticity at (i? = 3.1 a±,Z = 0). The densities have been 
normalized to one. 
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FIG. 6. Current field around the vortex core corresponding to the vortex ring state calculated 
taking = 0.2 uj± in Eq. (|l2|). The current is in arbitrary units and the axes are in units of a±. 
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